Low energy dynamics of the one dimensional multichannel Kondo-Heisenberg Lattice 
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We determine exactly the fixed point Hamiltonian of the one dimensional multichannel Kondo- 
Heisenberg lattice model for any number of channels N > 2. An anomalous singlet with non trivial 
internal dynamics is generated. We compute the correlation functions of the various conventional 
and unconventional order parameters of the system and find that for N < 4 the composite order 
parameter induce the dominant instabilities. 



The Kondo lattice is one of the most challenging 
problems in contemporary theoretical condensed matter 
physics. If the single impurity Kondo problem, a lo- 
cal moment antiferromagnetically coupled to conduction 
electrons, is by now well understood thanks to a vari- 
ety of theoretical techniques Jjj] the problem of a regular 
three-dimensional array of local moments in a metal (the 
Kondo lattice) still defies theoretical analysis. The diffi- 
culty stems from the fact that there are two competing 
effects, the tendency of the local moments to form sin- 
glets with the conduction electrons (or more complicated 
states H if more than one band is available) and the 
tendency of the moments to order due to the RKKY in- 
teraction mediated by the conduction electrons. Besides 
its intrinsic theoretical interest for the general theory of 
strongly correlated fermions, the Kondo lattice model is 
believed to capture the non- Fermi liquid physics of a class 
of rare-earth or actinide compounds |3| in which f-shell lo- 
cal moments couple to the electron bands. In this context 
H , it has been suggested that two-channel Kondo model 
could be relevant in Ce 3+ or U 4+ based alloys such as 
UBei3 and CeCu2Si2- This hypothesis is supported by 
experiments on diluted alloys that show a two-channel 
single impurity behavior in the specific heat or magnetic 
susceptibility. A full understanding of the multichannel 
Kondo lattice would also be of great significance for ex- 
periment. 

Many problems of strongly correlated fermions can be 
treated very effectively in one dimension thanks to pow- 
erful analytical [0 and numerical techniques available 
there. In particular, the full phase diagram of the one- 
dimensional single channel Kondo Lattice has been deter- 
mined Q] . An insulating state with a spin gap obtains at 
half filling. Away from half filling a small Kondo coupling 
produces a paramagnetic metal and a large one leads to 
ferromagnetism. Adding an antiferromagnetic Heisen- 
berg interaction between the local moments (Kondo- 
Heisenberg lattice) produces a spin-gapped metal Q with 
non-conventional pairing fluctuations Q . 

In this Letter, we investigate the ground state, excita- 
tions and dominant instabilities of the multichannel one 
dimensional Kondo-Heisenberg lattice model for incom- 
mensurate fillings and zero temperature. We find that 



screening of the local moments occurs via the " chiral sta- 
bilization" mechanism p"2[ , resulting in a ground state 
that is a chiral non-Fermi liquid we call coset singlet. 
We show that the dominant instability for a low num- 
ber of channels N < 4 is of the non-conventional pairing 
type as in the single channel case (|] and that for a larger 
number of channels, one recovers conventional pairing in- 
stabilities as the dominant ones. Fermi Liquid physics is 
regained in the limit of very large number of channels. 
We also discuss the effect of channel anisotropy. 

The Hamiltonian of the one dimensional multichannel 
Kondo-Heisenberg lattice is: 

H = —t (cl n <J Ci+i t n ya + cl +l n ty Ci y n,tr) 

(1) 

where a =t, I is the spin of an electron, n = 1, . . . , N 
is the channel index and Si is a localized spin 1/2. £ is 
the bandwidth, Xh > the Heisenberg coupling of the 
localized spins and Xk > the Kondo coupling. 

To study the low energy physics of the model we fol- 
low the standard strategy : keep only linear electron 
modes around ±kp captured in terms of the left and 
right moving fields V'l.-R.ct.^x), as well as the low-lying 
spin-lattice modes around -. The low energy physics is 
then described by a continuum Hamiltonian expressed 
using non-abelian bosonization [jl0|ll| as a quadratic 
form of various currents: the SU{2)\ (left and right) 
spin currents <jl,r of the local moments, the SU(2)n 
(left and right) spin current Sl,r(x) of the electrons and 
the electron charge and channel currents which decouple 
from the spin currents in the continuum. The notation 
SU(M)k indicates that the currents (generically denoted 
J1 R ) satisfy a Kac-Moody algebra: [J£(x), </£(x')] = 
5{x - x')f abc Jl(x) + i<5'(x - x'), where f abc are the 
structure constants of the SU(M) Lie algebra. A similar 
relation holds for the currents Jr, while the Jr and Jt, 
currents commute. The KM algebra generates a confor- 
mal field theory (CFT) known in Lagrangian form as the 
Wess-Zumino-Novikov-Witten (WZNW) model. All 
the physical operators of the original lattice theory (|l|) 
can be expressed in terms of the primary or descendant 
fields of this CFT allowing the calculation of the asymp- 
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totic behavior of their correlators. This is at the heart of 
all applications of non-abelian bosonization to condensed 
matter physics. 

Consider the Hamiltonian with Xk — 0. The elec- 
trons are free, and in the basis described above charge, 
spin and channel excitations propagate independently of 
each other, and decouple from the local-moment exci- 
tations. Turning on a small antiferromagnetic Kondo 
coupling < Xk -C t, Xh at incommensurate filling, the 
charge and flavor excitations remain decoupled from the 
spin excitations. We can thus restrict ourselves to the 
latter, described by the Hamiltonian: 
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(Sl-Sl + Sr.Sr) H ^-{<Jl-&l + vr-vr) 
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+ X f K (S L .a L + S R .a R ) + X b K (S L .a R + S R .a L )] (2) 



where we dropped irrelevant terms such as <jl-3r The 
coupling X K describes the forward scattering, X h K the 
backward scattering and X K — X b K — Xk to begin with. 
Under RG transformations X K does not flow near the 
weak coupling fixed point and merely renormalizes the 
velocity. On the other hand X h K is relevant and drives 
the system to a new strong coupling fixed point. We will 
thus take X K = and v s = v in (||) and determine the 
fixed point. We will later prove that the neglected terms 
are irrelevant near the strong coupling fixed point. 

Under these circumstances the spin sector is described 
by H = Hi + H2 where: 



2nv ? - 2ttu , x ? - 

N + 2 — i~ aL<JL * K ^ R - aL 

2™ g =? 2nv ^ _ ~ _ ' 

N + 2 — ^-vrVr + Xk^l-^r. 



Hi= J dx 
H 2 = f dx 



In Hi (resp. H2), the left (resp. right) branch of the 
SU{2) N WZNW model is coupled to the right (resp. left) 
branch of the SU{2)\ WZNW model, leading to a chiral 
asymmetry of Hi (resp. H 2 ). We readily identify the 
strong coupling fixed point of (||) exactly via "chiral sta- 
bilization" [ ^2[ . The chiral asymmetry in Hi or H 2 is 
invariant under the RG flow and characterizes the fixed 
point. We find this way that the theory (|J) flows under 
RG to a fixed point theory that is the product of a coset 
theory Q by a WZNW theory: 

bu (2) N 

where the coset theory, s,7 ( 2 )^^^ 2 )w-i describes the 
local moment spin sector, and the SU(2)n-i WZNW 
theory describes the electron spin sector. Note that at 
the fixed point the left and right components is H\ and 
H2 are recombined and chiral symmetry is globally pre- 
served. 



What is the physics around the fixed point? The coset 
theory describes a spin singlet which the local moments 
form with the electrons. It is a new type of a singlet, a 
coset singlet: a fraction ( N+ i^ N+2 ) of the local moment 
"modes" are paired with the same number of electron 
spin "modes" . Thus the system loses twice this amount 
of degrees of freedom as seen in the total specific heat 
(including channel and charge degrees of freedom): 



Ctotal - + 1 - {N + ^ 2 iN + 2) 

The susceptibility is given by, 



T. 



(4) 



(5) 



2N+1 12 

and the Wilson Ratio: R w = MEMM , The 

coset singlet still retains degrees of freedom whose num- 
ber is given by the central charge of the theory c = 



1 



This fraction decreases with the number 



(N+l){N+2)- 

of channels N since it is "easier" to form the singlet when 
N increases. This also shows up in the effective coupling 
of the electrons to the local moments which decreases 
with with the number of channels, At- ~ 1/N. 

Consider the two-channel case, ]14| . For N = 2, lo- 
cal moments are described by a ^^g^^^ 1 = Ising 
theory, or equivalently by a Majorana fermion. Such a 
Majorana fermion picture is very appealing since it is 
well known pTEj j that the the single impurity two chan- 
nel Kondo model be described at the fixed point by a 
local Majorana fermion degree of freedom. These Majo- 
rana fermions form a band when coupled with each other, 
thus suppressing the single impurity residual entropy at 
T = 0. 

Having obtained the low energy theory ([|) describing 
the spin sector, we proceed to express the original oper- 
ators in terms of the operators of the fixed point theory. 
This will enable us to check that the operators we dis- 
carded in (|J) are indeed irrelevant at the strong coupling 
fixed point as well as determine the dominant instabil- 
ity of the Multichannel Kondo-Heisenberg Lattice. The 
needed identification of operators as well as the calcula- 
tion of the scaling dimensions was done in Ref. [O] . Let 
us summarize briefly the method. To obtain the con- 
formal weight of a given operator, we first decompose it 
into a product of operators belonging to each of the two 
decoupled chiral theories. Then, for each chiral theory, 
we decompose operators of SU(2)n ® SU(2)i on oper- 
ators of SU(2)n_i in an expansion formally similar to 
the Clebsch-Gordan expansion, the role of the Clebsch- 
Gordan coefficients being played by operators of the coset 
theory jl2| . The operator with the lowest scaling dimen- 
sion in this expansion is then retained as the fixed point 
form of the original operator. 

The results are summarized in the table | for the theory 
described by H±. The conformal weights of the theory 
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described by H 2 are obtained by interchanging L and 
R. These conformal weights are such that the operators 
S L (x).a L (x), S R (x).a R (x) and (v s - v)(a R .a R + a L .a L ) 
that we have previously discarded are indeed irrelevant 
(marginally irrelevant for N = 2). This proves the self- 
consistency of our treatment. 

The fixed point is a non-Fermi liquid. The Green's 
functions of the right moving fermions is given by: 



(T^(x,t)^(0,0)> 



1 



(x - vt) i+ 2 s "( x + vty. 



Here we have taken, for sim- 



where S N = {N+1){N+2) . 
plicity, charge, spin, and channel velocities to be equal, 
and combined exponents from all sectors, neglecting the 
non-universal Luttinger interaction in the charge sector 
(it can be easily taken into account JT2|| ). Further, if a 
contribution in the channel sector is generated it is fer- 
romagnetic and flows to zero. GL(x,t) is obtained by 
replacing x ± vt by x =F vt. We note that a weak sin- 
gularity appears at the Fermi level Uf 1 and there is no 
large Fermi surface. Also note that all dimensions tend 
to their Fermi liquid values in the limit of large number 
of channels. 

We now examine the possible order parameters of 
the system. Beginning with the localized moments, 
we have Si — <tl(x) + a R {x) + e l7TX ^ a n(x), where 
n{x), the staggered magnetization, is given by n{x) — 
\ 9 r a ^a.(39L,i3, where g Rta , (g L . a ) are the right 
(left) WZNW fields. We find, 



(n(x).nV)) 



1 



(G) 



More order parameters are available in the electron 
sector. The order parameters for charge density wave 
(CDW), spin density wave (SDW), singlet (SS) and 
triplet (TS) superconducting are defined in the Table and 
their dimension is given, from which it follows that they 
are degenerate and fall with a power 2 + 2<5jv- As the 
fluctuations of these order parameters are weaker than 
in the one dimensional metal we are lead to investigate 
the possibility of dominant fluctuations associated with 
a non-conventional order parameter j|] odd-frequency 
singlet pairing [|l6| (c-SP) and composite charge-density 
wave order (c-CDW): O c -SP = n.Ors an d O c -cdw = 
"•Osdw- The composite operators have momentum ±^ 
and 2kp±^ respectively, and are associated with the gap- 
less excitations predicted by Yamanaka et al. |17|j ; their 
correlations are expected display quasi-long range order. 
Indeed, they decay with the power 3 — j^p^- 

We observe that for N < 5, the composite order pa- 
rameters have the most divergent correlations. In this 
case a large enough fraction of electron spin degrees of 
freedom is bound the the local moments to suppress the 
conventional order parameters and enhance the compos- 



ite ones. This situation is similar to one-channel one di- 
mensional Kondo Lattice case || where composite pair- 
ing operators are also dominant. For N — 5, the two 
types of order are degenerate at the fixed point, and one 
needs to describe also the approach to the fixed point 
starting from the bare Hamiltonian to determine the or- 
der. For N > 6, the fraction of electron spin degrees 
of freedom bound to a local moment is insufficient to 
permit composite order parameters to dominate the con- 
ventional ones, thus recovering the Fermi liquid limit. 

Thus far, we assumed channel isotropy. In the sin- 
gle impurity problem Q channel anisotropy is a relevant 
perturbation, so we have to investigate its effect for a lat- 
tice. Assume Ni channels couple to the local moments 
with coupling strength A^ and N 2 — N — Ni channels 
couple with strength X 2 K . Then, the spin excitations of 
the N\ channels are described by a SU{2)m 1 KM alge- 
bra, whereas those of the N 2 remaining channels are de- 
scribed by a SU(2)n 2 . For <C \ 2 K , coset screening 
occurs between the N\ channels and the local moments. 
The resulting theory is: ^g^^" 1 ' 1 ® SU(2) Nl ^ x 
SU(2) m 21 with central charge: c = 2N + 1 — (jVi+i)(iVi+2) 
smaller than c = 2N + 1 — (jv+i)(JV+2) ^ e c h ar S e °f the 
symmetric fixed point. It is the stable fixed point having 
the lower number of degrees of freedom []l8|| . This implies 
that channel anisotropy is always relevant as in the single 
impurity problem. 

We have derived the fixed point theory describing the 
one dimensional multichannel Kondo-Heisenberg lattice. 
This fixed point is in the class of chiral non-fermi liq- 
uids. We showed that the dominant instability is of the 
composite pairing type for a number of channels smaller 
than five and of the conventional pairing type otherwise. 
Our results were obtained in the limit of weak Kondo 
coupling, at zero temperature and away from half-filling. 
At half-filling, it is known that a spin gap develops only 
for N-2S integer 

There are many directions for future research: At fi- 
nite temperature, the behavior we described must disap- 
pear above the larger of the Kondo Temperatures T^- 
and Tjf characterizing the electron and local moment 
sectors, respectively. Also the problem of exhaustion 
pl[ , how a large concentration of impurities reduces the 
Kondo scale, is still unresolved. The conformal field the- 
ory approach of our paper is not able to settle this issue. 
However, the model (||) happens to be Bethe Ansatz in- 
tegrable. This will allow us to settle the issue of the 
Kondo scales and their dependence on the filling - ex- 
haustion - and discuss the full crossover to the zero tem- 
perature behavior. A simple generalization is to replace 
the spin-1/2 chain by an integrable spin-S* chain of lo- 
cal moments, which would lead to a replacement of the 
SU(2)i by SU{2) 2 s Kac-Moody algebra. Another line 
of future research relates to coupling two or more Kondo 
chains to investigate the effects of the unstable fluctu- 
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ations. A last question is whether the picture we have 
obtained persists for \k 3> t, \h or very low density and 
if not what is the strong Xk regime. Presumably, such 
regime would be a ferromagnetic Nagaoka-like state as 
in the single channel case. Besides the various general- 
izations of the one dimensional Kondo-Heisenberg lattice 
problem, it would be interesting to determine whether 
the physics of the Kondo-Heisenberg problem persists in 
the Kondo limit Xh = 0. It is known that this is not 
so in the one-channel case 0,9- However, since in the 
multichannel case there is only a partial screening of the 
electrons by the spins, one may expect that a RKKY 
interaction could be generated even in the pure Kondo 
problem putting it in the universality class of the Kondo- 
Heisenberg problem. 
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TABLE I. The conformal weights of the operators in the 
theory described by (^|). Here, Sn = (jv+i) 3 (iv+2) • F° r N — 2, 
the conformal weights of Sl and or are respectively (0, 1) and 



(1,0). 




Operator 


Conformal weights 
at the fixed point 

11 i 6w d N \ 



9L,0 
9R,0 

a L (x) 
or(x) 
Sr(x) 
S L (x) 

^ OcDW = i>t, n ,^R,n,a 
OsDW = W i L, n ,<j°a,a'' l i'R,n,v 
Oss = -i1pt,n,affl i<T ,ll>R,n,o 
Ots = -~Ml>L,n, a (pa y ) a ^ l ^>R. 

Oc-SP = U.OtS 

Oc-cdw = n.OsDw 



2 ' 4 ' 4 

4 > 2 ' 4 

3 1 i 3 



4(JV+1) ' 4 



4(iV+l) 
3 



4(JV+1) ' 4(JV+1) 



I AT+1 ' 1 d" jv+l 

(1 + -2 2- 

V ' N+l ' AT+1 , 



,N > 3 
, N > 3 




2(JV+2) 
3 



2(JV+2) ' 4 2(iV+2) 
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